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DEVELOPMENT OF SKYLAB EXPERIMENT T-013 
CREW/VEHICLE DISTURBANCES 


By Bruce A. Conway 
Langley Research Center 

SUMMARY 

A Skylab experiment to determine the characteristics and effects of crew-motion 
disturbances has been developed. The experiment, designated T-013, will correlate data 
from histories of specified astronaut body motions, the disturbance forces and torques 
produced by these motions, and the resultant spacecraft control-system response to the 
disturbances. Primary application of crew-motion-disturbance data will be to the sizing 
and design of future manned spacecraft control and stabilization systems. This report 
presents the development of the crew/vehicle disturbances experiment. Appendixes 
derive a mathematical model of human body motion which may be used for analysis of a 
variety of man-motion activities. 


INTRODUCTION 

Motions of astronaut crew members while in contact with a spacecraft can produce 
disturbances on the spacecraft which would affect pointing accuracy, or otherwise impose 
demands on the stabilization and control system. The resultant attitude errors (ref. 1) 
can be cumulative if the crew -motion path results in a net rotation of crew mass relative 
to the spacecraft center of mass. In addition, any transient rate and attitude errors are 
generally sufficiently large to require active control. Sizing of such a control system 
is dependent on the crew-motion disturbance characteristics. Since crew motions gen- 
erally represent the largest disturbance torques and forces for manned spacecraft, and 
the frequencies contained in these torque histories are several orders of magnitude 
higher than those associated with external disturbances (such as gravity -gradient and 
aerodynamic moments), the torque level and bandwidth of the control system must be 
sufficiently large to allow compensation for the crew disturbances. Flexibility modes 
of the spacecraft and any attached experiment packages may also be excited by crew 
motions, necessitating definition of frequency content of the disturbances. Simple extrap- 
olation of the sizing requirements of the control system for one-, two-, or three -man 
crews to account for six-, 12 -, or even 50-man crews might result in unnecessarily large, 
overdesigned systems. Detailed characteristics of crew-motion disturbances including 
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both discrete and stochastic models of crew-motion forces and torques must thus be 
developed to permit effective design of future spacecraft and experiment control systems. 
Accurate knowledge of these characteristics will result in realistic control-system sizing 
requirements. 

Initial studies of crew-motion disturbances considered crew members as point 
masses, and investigated the effects on a spacecraft of translating one or more of these 
masses from one point to another within the spacecraft. The studies included digital 
computer simulations of spacecraft dynamics, along with spin tests using dynamically 
scaled models of various spacecraft inertia configurations. Following these earlier 
simulations, consideration of man, himself, as a dynamic system was begun. Develop- 
ment of an anthropometric model of the human body was carried out, using earlier U.S. 

Air Force anthropometric studies (refs. 2 to 6). This anthropometric model was extended 
to include the dynamics of body-segment motions. The resultant analytical model may be 
used to compute crew-disturbance forces and torques in a spacecraft and the spacecraft 
response, for specified astronaut motions. Efforts carried out by Hendricks, Johnson, 
and Kirch (ref. 7) have developed stochastic representations of crew -motions -disturbance 
forces and moments for selected crew activities, with applications to control-system 
frequency-response investigations. 

Use of both discrete and stochastic techniques developed to describe crew-motion 
disturbances requires verification of these techniques by experimentation. Several tech- 
niques of simulating the zero-gravity environment of a manned spacecraft and their pri- 
mary advantages with respect to a crew-motion experiment are listed in table I. Because 
of the limitations of ground-based simulations and, consequently, the questionable appli- 
cability of resultant data, a manned space flight experiment was proposed to be flown 
aboard the Skylab spacecraft. The Skylab orbital workshop (see fig. 1) provides an excel- 
lent laboratory for crew disturbance research, since it is primarily an experiment car- 
rier, and hence has facilities to provide maximum experiment support. 

This report will discuss the development of a manned spaceflight experiment whose 
objectives are (1) to determine the characteristics of typical crew-motion disturbances 
and (2) to assess the response of a representative manned spacecraft control system to 
some of these disturbances. The resultant experiment, designated Skylab experiment 
T-013, will be described with respect to its implementation hardware and protocol. 

Included in the report as pertinent to the development of experiment T-013 are two 
appendixes: Appendix A, by Charles T. Woolley, presents the derivation of the physical 
properties for the anthropometric model of man; appendix B, by Peter R. Kurzhals and 
Robert B. Reynolds, derives the dynamic equations associated with motions of the anthro- 
pometric model. 
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TABLE I.- SEVERAL ZERO-GRAVITY SIMULATION TECHNIQUES 


Simulation approach 
Underwater neutral buoyancy 

Advantages 

Disadvantages 

Approximates actual zero-g 
for unsuited subject 

Drag excessive for all 
but slowest motions 
Breathing equipment 
restrictive 

Air bearing simulation 

i 

Provides good low- or 
zero-g effect in two 
dimensions 

Implementation of force 
measuring techniques 
difficult 

Requires extremely fine 
balance and c.g. shift 
compensation 
Must counteract gravity in 
many motions 
Mounting harnesses, etc., 
too restrictive 
Susceptible to ambient air 
movement 

Servo-drive simulation 
(ref. 8) 

Can be tied together 
directly with computer 
simulation of spacecraft 
dynamics 

Must counteract gravity in 
many motions 
Mounting harnesses too 
restrictive 

Cable suspension 

Relatively low cost 

Degrees of freedom limited 
Pendulum effects present 
Support apparatus 
restrictive 

Zero-g aircraft 

Actual zero-g environment 

Short run times 
Unnatural positive -g forces 
interspersed between 
zero-g runs 
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Euler angle transformation matrix 
element of [d] 

lower torso hinge point location with respect to bottom of torso 
of model man 

force 

,Fg force in load cell 1,2,. . .,6 
angular momentum 
moment of inertia 

location of load-cell-array centroid 

moment 

mass 

corrected total mass of model man 
mass factor (see eq. (Bll)) 
position radius vectors (see fig. 36) 


RR 

SL 

SR 

SML 

SMR^ 




model -man segment dimensions defined in appendix A 
(see figs. 18 to 30) 


SIxX>SIyy,SI Z z model-man segment moments of inertia 


S.D. standard deviation 


T 


torque 



t 


time 


V volume 

X,Y,Z coordinate axes of spacecraft body 
Xp,Yp,Zp inertially fixed axes 

Xj,Yj,Zj intermediate spacecraft axes (see fig. 36) 

X M ,Y M ,Z M model-man center -of -mass body axes 

X MI’ Y MI’ Z MI model -man reference axes 

x MS’ Y MS> z MS model-man center -of -mass axes parallel to spacecraft body axes 

x,y,z Cartesian coordinates 

y model-man segment Euler angle 

V fraction of frustum height at which frustum center of mass is located 

|U ratio of base radii of conical frustum 

77 = 3.14159 

p density, mass per unit volume 

cr = 1 + jii + 

^,8,cp Euler angles for spacecraft (see fig. 39) 

fl angular velocity 

spacecraft body rates 

Subscripts: 

cal related to calibration values 
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related to center of mass 

external 

foot 

head 

hand 

jth segment of model man (j = 1, 2, . . 

lower arm 

lower leg 

lower torso 

man 

related to reference axes of model man 

related to geometric spacecraft origin 

rigid body 

spacecraft 

total 

upper arm 
upper leg 
upper torso 

with respect to X-, Y-, and Z-axis 
segment 1,2,. . 9 of model man 



Dots over symbols indicate derivatives with respect to time. 


O 
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denotes a vector quantity (see appendix B) 
denotes a matrix (see appendix B) 

EXPERIMENT HARDWARE DEVELOPMENT 


Determination of the characteristics of crew-motion disturbances, including the 
verification of analytical crew-motion models, requires that certain data be obtained 
from experiment T-013. Body-limb motion and center-of-mass translation histories 
must be measured simultaneously with the measurement of any forces and moments pro- 
duced by the motions. During portions of the experiment, vehicle control-system 
response must also be measured simultaneously in order to assess its performance in 
the presence of known crew disturbances. The following sections will discuss the defini- 
tion and development of flight hardware to meet the objectives of the Crew/Vehicle 
Disturbances Experiment. 


Limb-Motion Sensors 

The two prime candidates for a limb-motion sensor system (hereinafter referred 
to as LIMS) evaluated during the definition of experiment T-013 were motion pictures and 
an exoskeletal device. Trade-off studies were performed to determine whether an exo- 
skeleton or a photographic system would be more advantageous in the flight experiment 
limb-motion measurement. The prime advantage of motion pictures is that they permit 
complete freedom of body motion, although certain restrictions might be imposed on the 
experiment to ensure that all necessary limb segments would be in the camera's field of 
view. In addition, motion pictures could utilize existing onboard cameras for the limb- 
motion measurement function, as well as providing astronaut center-of-mass location 
and torso attitude data. Their outstanding disadvantages are the large amount (980 m) 
of film required to obtain necessary data resolution and the extensive amount of manual 
data reduction needed to obtain the time histories of the body-segment motions. Advan- 
tages of an exoskeleton LIMS, utilizing potentiometers at principal rotation points on the 
human body, are the ease of data reduction and the resultant savings in experiment weight. 
Its major disadvantage is that some restriction is imposed on subject motion; however, 
proper design and sizing of the exoskeleton could minimize its restrictions so that virtu- 
ally all routine body movements could be performed in a natural manner. Table II sum- 
marizes the results of the trade-off comparison. It was concluded that the exoskeleton 
offers advantages over motion pictures and this concept was selected. 
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TABLE n.- TRADE OFF BETWEEN EXOSKELETAL AND MOTION -PICTURE 
CONCEPTS FOR LIMB MOTION MEASUREMENTS 



Exoskeleton 

Motion pictures 

Orbited weight, kg 

*4.2 

6.4 

Orbited volume, 

*0.11 

0.02 

Data reduction cost per minute of 

experiment time (approximate) 

<$20 

>$800 

Return weight, kg 

0 

6.4 

Return volume, m^ 

0 

0.02 

Restricts body motion 

Possible 

Could impose 

Measurement accuracy, 

degrees of limb motion 

±1.2 

restraints for 
camera coverage 

±3.2 


A laboratory model LIMS (shown in fig. 2) was constructed and used to verify its 
measurement compatibility with the dynamic equations given in appendix B. Optimized 
design concepts for a lightweight exoskeleton LIMS were developed. These concepts 
included the mating of the exoskeleton to a one-piece form -fitting garment (similar to 
coveralls) that would provide for the maintenance of alinement between exoskeleton and 
body limbs and would simplify donning and doffing. This one-piece suit would be of such 
material and construction that a stretch capability would exist in order to allow for 
possible variations in experiment subject sizes. Figure 3 depicts the LIMS to be used 
in the flight experiment. It consists of an aluminum exoskeleton, with 16 potentiometers 
to measure rotations at the principal body joints (see fig. 4). Power (at 5 V d.c.) to the 
potentiometers and data signal return is via an electrical umbilical attached to the LIMS 
and the experiment data system. The exoskeleton is adjustable and can fit the entire 
range of astronaut sizes, while expansion panels in the one-piece suit extend its range of 
fit to several subject sizes. The LIMS, along with its data umbilical, is launched in a 
stowage container (shown in fig. 5). This container provides protection from the loads 
imposed during launch of the spacecraft and also affords protection from dust and debris 
during periods of nonuse. 


Force Measuring System 

The measurement of forces and moments produced by the T-013 subject astronaut 
requires apparatus similar to the six-component strain-gage balances used in many wind 
tunnels. However, it must be rugged enough to withstand normal handling and launch 
vibration loads. An array of load cells was selected as being the easiest to fabricate and 
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maintain (with respect to a strain-gage balance, for example), as well as having the 
capability to survive the launch environment. Of the load-cell array configurations con- 
sidered, the configuration shown schematically in figure 6 was selected. In addition to 
providing minimum weight and volume, this configuration provides essentially equal sen- 
sitivity about all three measurement axes. Using this arrangement, triaxial forces and 
moments are computed from the following relations: 

F x = cos 45°[F 1 - F 6 + cos 60°(F 2 - Fg + F 4 - F 5 )J + F xcal 



Figure 4 . - LIMS exoskeleton showing potentiometer locations. 
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BASE " 

Figure 6.- Schematic drawing of load-cell configuration selected 
for flight hardware implementation. 

other and to provide a measure of redundancy in that either could he used for evaluation 
of activity where the subject is attached to a single location in the spacecraft. Foot 
restraints (shown in fig. 8) on one FMU permit the astronaut to be restrained for studies 
of "restrained" or "in-place" motion effects. 

Astronaut Position and Attitude Determination 

In order to define completely an astronaut’s location and orientation within a space- 
craft, histories of his center-of-mass motion and torso attitude must be known along with 
his limb motions, since all limb-motion measurements are referenced to the torso (or 
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Figure 7.- Laboratory model of foi 


input and data readout e 




L-71-7125 


measuring unit. Associated power 
:tronics are shown also. 





segment 1 of the analytical model). Several means of monitoring these histories during 
the performance of experiment T-013 were investigated. These methods include slow- 
rate motion pictures, inertial measuring unit, and ultrasonic audio signal tracking. 

The motion-picture concept utilized two cameras and conventional photographic 
triangulation. Frame rate of the cameras was 6 frames per second, which would be 
adequate to obtain the necessary torso attitude and center-of-mass translation histories, 
since expected rates of motion for these quantities is significantly lower than rates of 
individual limb articulations. 

The inertial-measuring -unit concept consisted of six accelerometers mounted on 
an orthogonal yoke worn by the astronaut. This scheme would require periodic updating 
(whenever acceleration levels dropped below sensor threshold, for example) and would 
need highly accurate transducers to be able to determine distance with the same accuracy 
as velocity. 
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The ultrasonic -triangulation method used an ultrasonic transmitter initiating a 
pulse simultaneously with a radio frequency (RF) transmitter, both located on the subject. 
The delay time between receiving the RF pulse and the ultrasonic pulse would be directly 
dependent on the distance from the receiver to the astronaut. This method would require 
multiple receivers and several transmitters in order to determine both center-of-mass 
location and body torso orientation. It would also require knowledge of the spacecraft 
atmospheric properties for determination of the speed of sound in that medium. 

The photographic technique was chosen because it made use of existing onboard 
camera equipment, and did not entail the mechanization complexity of the other methods. 
Appendix G of reference 8 presents an example of a method of photogrammetry that may 
be used to reduce T-013 astronaut position and attitude data. 


Experiment Data System 

The experiment data system is designed in such a manner that data reduction is 
facilitated. For T-013, measurements must include data from the LIMS and FMU's, 
along with a timing reference from the spacecraft central timing system. A commutator, 
which would sequentially sample each experiment transducer and the timing signal, was 
selected. The commutated data will be formatted into a single serial bit stream for input 
to the spacecraft tape recorders and stored for subsequent telemetry to ground receiving 
stations. Figure 9 shows the flight data system which features a switch to activate power 
to the LIMS, through the data umbilical. Table III lists the experiment measurements to 
be processed by the T-013 data system. 
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TABLE HI.- MEASUREMENTS TO BE PROCESSED BY T-013 DATA SYSTEM 


Skylab 

measurement 

number 

Measurement name 

Sample rate, 
per sec 

G7017 T-013 

LIMS angle, right shoulder, rear 

20 

G7023 T-013 

LIMS angle, right shoulder, forward 

20 

G7024 T-013 

LIMS angle, right upper arm 

20 

G7055 T-013 

LIMS angle, left shoulder, rear 

20 

G7054 T-013 

LIMS angle, left shoulder, forward 

20 

G7021 T-013 

LIMS angle, left upper arm 

20 

G7022 T-013 

LIMS angle, right elbow 

20 

G7057 T-013 

LIMS angle, left elbow 

20 

G7058 T-013 

LIMS angle, right hip, rear 

20 

G7056 T-013 

LIMS angle, right hip, forward 

20 

G7059 T-013 

LIMS angle, right thigh 

20 

G7060 T-013 

LIMS angle, left hip, rear 

20 

G7061 T-013 

LIMS angle, left hip forward 

20 

G7063 T-013 

LIMS angle, left thigh 

20 

G7062 T-013 

LIMS angle, right knee 

20 

G7052 T-013 

LIMS angle, left knee 

20 

S7005 T-013 

FMU #1, load cell 1 

20 

S7004 T-013 

FMU #1, load cell 2 

20 

S7003 T-013 

FMU #1, load cell 3 

20 

S7007 T-013 

FMU #1, load cell 4 

20 

S7006 T-013 

FMU #1, load cell 5 

20 

S7009 T-013 

FMU #1 , load cell 6 

20 

S7008 T-013 

FMU #2, load cell 1 

20 

S7001 T-013 

FMU #2, load cell 2 

20 

S7000 T-013 

FMU #2, load cell 3 

20 

S7011 T-013 

FMU #2, load cell 4 

20 

S7010 T-013 

FMU #2, load cell 5 

20 

S7002 T-013 

FMU #2, load cell 6 

20 

M7075 T-013 

FMU #1, calibration voltage 

10 

M7076 T-013 

FMU #2, calibration voltage 

10 

M7077 T-013 

LIMS, calibration voltage 

10 

K7328 T-013 

Airlock module timing reference 

10 

K7329 T-013 

T-013 identification word 
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Spacecraft Pointing Control System 

An important adjunct to the Crew/Vehicle Disturbances Experiment is the evalua- 
tion of spacecraft-control-system response to known crew -motion -disturbance inputs. 
Data from the Apollo telescope mount (ATM) pointing control system (described in ref. 9) 
will be obtained during the performance of experiment T-013 and will include control- 
moment gyro (CMG) gimbal angle and rate measurements, rate gyro measurements, sun 
sensor signals, and control-moment gyro command parameters. A complete list of ATM 
measurements to be obtained is found in table IV. This data will be subsequently time- 


TABLE IV.- ATM MEASUREMENTS TO BE OBTAINED DURING EXPERIMENT T-013 


Skylab 

measurement 

number 

Measurement name 

Sample rate, 
per sec 

K122 

ATM mission elapsed time 

10 

R017-702 

CMG #1, outer gimbal rate 

12 

R018-702 

CMG #2, outer gimbal rate 

12 

R019-702 

CMG #3, outer gimbal rate 

12 

R020-702 

CMG #1, inner gimbal rate 

12 

R021-702 

CMG #2, inner gimbal rate 

12 

R022-702 

CMG #3, inner gimbal rate 

12 

M164-702 

Fine sun sensor, yaw output voltage 

12 

M165-702 

Fine sun sensor, pitch output voltage 

12 

R023-702 

Vehicle rate gyro No. 1, X-axis output 

12 

R024-702 

Vehicle rate gyro No. 1, Y-axis output 

12 

R025-702 

Vehicle rate gyro No. 1, Z-axis output 

12 

R026-702 

Vehicle rate gyro No. 2, X-axis output 

12 

R027-702 

Vehicle rate gyro No. 2, Y-axis output 

12 

R028-702 

Vehicle rate gyro No. 2, Z-axis output 

12 

R029-702 

Vehicle rate gyro No. 3, X-axis output 

12 

R030-702 

Vehicle rate gyro No. 3, Y-axis output 

12 

R031-702 

Vehicle rate gyro No. 3, Z-axis output 

12 

G023-702 

CMG #1, inner gimbal angle 

12 

G024-702 

CMG #1 , outer gimbal angle 

12 

G025-702 

CMG #2, inner gimbal angle 

12 

G026-702 

CMG #2 , outer gimbal angle 

12 

G027-702 

CMG #3, inner gimbal angle 

12 

G028-702 

CMG #3, outer gimbal angle 

12 

R015-702 

Experiment pointing system pitch rate 
gyro command 

12 

R016-702 

Experiment pointing system yaw rate 
gyro command 

12 

M529-702 

Acquisition sun sensor, X output voltage, 
primary 

12 

M530-702 

Acquisition sun sensor, X output voltage, 
redundant 

12 

M531-702 

Acquisition sun sensor, Y output voltage, 
primary 

12 

M532-702 

Acquisition sun sensor, Y output voltage, 
redundant 

12 

G037-702 

Star tracker, outer gimbal angle 

12 

G038-702 

Star tracker, inner gimbal angle 

12 

K382-702 

Experiment pointing control system resolver 
output (pitch and yaw) 

1.25 
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correlated with the FMU and LIMS data. The measurements listed in table IV will be 
telemetered in real-time while the orbital workshop (OWS) is over ground receiving sta- 
tions, since available ATM data recorders record at a rate too low for adequate data 
resolution, and do not have the capability to record all measurements listed in table IV. 

EXPERIMENT DESCRIPTION 


Experiment T-013 will be performed in the forward area of the Skylab orbital work- 
shop. Figure 10 depicts the layout of this section of the spacecraft, and shows the loca- 
tion of the various pieces of T-013 experiment hardware (including camera placements). 



CAMERA (OWS DOME) 



Figure 10.- Sketch of Skylab spacecraft interior showing locations of 

T-013 experiment hardware. 

As noted on this figure, the two force measuring units are approximately 2.75 meters 
apart, with their sensing plates parallel, and are attached to the OWS wall structure via 
permanently installed brackets. 

The orbital assembly (which includes the orbital workshop and the docked Apollo 
command-service module) will be in a nominal 435-km circular orbit, with an orbit 
inclination of 50°. The orbital assembly will be attitude -stabilized in the solar inertial 
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mode during T-013. This mode provides for keeping the Apollo telescope mount pointed 
toward the sun, with the long axis of the orbital workshop maintained in the orbit plane. 
The control moment gyros of the ATM pointing control system (PCS) will provide the 
primary means for spacecraft stabilization. 

The T-013 subject will perform a variety of motions during the Crew/Vehicle 
Disturbances Experiment. These will be restricted as far as possible to motions of the 
9-segment body described in the appendixes. One set of motions will be conducted while 
attached to one EMU; these motions include simple arm and leg pendulum -type move- 
ments, random reaching and twisting movements, and exaggerated breathing, coughing 
and ’’sneezing 1 ’ functions. Figure 11 shows an artist’s concept of this phase of the Skylab 
experiment. Free soaring translation between the FMU’s, with push-offs by either hands 
or feet, will comprise the remainder of the experiment activity. Certain atypical motions 
will be performed by the subject alone and, concurrently, with a second astronaut, in 
order to provide a "worst-case" crew -disturbance input to the Skylab pointing control 
system. These motions will consist of vigorous limb movements and rapid soaring 
translations. 







Most of the discrete motions will be performed three times by the subject, and the 
discrete -motion sequence will be repeated (for a total of six performances per motion) 
to give a degree of statistical confidence to the resultant data. The experiment perfor- 
mance timeline is summarized as figure 12. 

An observer astronaut will assist the subject in donning and doffing the LIMS and 
in attaching him to an FMU. This observer will enter voice comments into the tape 
recorder during experiment progress; these comments will include task identification 
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DEACTIVATION 

AND 

STOWAGE 


Figure 12.- T-013 experiment timeline summary. 

and relayed assessments by the subject as to difficulty in performing a given task or 
other operational peculiarities. Another responsibility of the observer will be to actuate 
the data cameras and tape recorder. 

Scheduling of the experiment during the Skylab flight will depend in part on avail- 
able ground tracking station coverage for the n worst-case n motion sequence and other 
portions of T-013. Typical telemetry coverage during a given station overpass is 
approximately 5 minutes. Figure 13 shows the ground tracks of three representative 
Skylab orbits that could be used for T-013. 
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Figure 13.- Skylab western hemisphere telemetry receiving station coverage, 
showing ground tracks of three representative Skylab orbits. 

CREW-MOTION DATA USE AND APPLICATIONS 

The data obtained from the Skylab Crew/Vehicle Disturbances Experiment will be 
used in the sizing, design, analysis, and simulation of control-system concepts and con- 
figurations for future manned space vehicles. Previous simulations of crew-motion dis- 
turbances have been conducted. One simulation (ref. 10) utilized a three -degree -of - 
freedom cable suspension scheme, with forces and moments measured by a strain-gage 
balance. Agreement between results of this simulation and predicted results using the 
analytical model of appendix B was poor. A second simulation (ref. 11) used a servo- 
drive six-degree -of -freedom simulator, and load-cell arrays similar to those for the 
T-013 flight hardware. Results of this simulation when compared to the mathematical 
model prediction (ref. 8) were in much closer agreement than the cable suspension results, 
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Figure 14. - Skylab experiment T-013 data use and applications. 



















although significant differences did exist. Variations in force and moment levels between 
experiment and theory are attributed to 1-g effects of the simulations, as pointed out in 
table I. It is expected that reduction of the data from T-013 will, in addition to providing 
verification of the analytical crew-motion model, permit the derivation of correction fac- 
tors which can be applied to 1-g simulations to approximate more closely zero-g crew- 
disturbance effects. The verified analytical model may be applied to investigation of 
man -motion effects, with respect to astronaut maneuvering units, work capability and task 
analysis, and optimized reorientation maneuvers of a free-floating astronaut, for example. 
Crew-disturbance data from T-013 will also provide a base for developing stochastic 
models of crew-motion disturbances, which may be applicable to the analysis and simu- 
lation of multicrew, long-term flight activities (as initiated in ref. 7, for example). Data 
from the ATM PCS will provide additional verification of the dynamic model derived in 
appendix B. These data will also aid in evaluating the effectiveness of the Skylab control 
system in compensating for crew disturbances and in providing high-accuracy pointing 
levels for a manned spacecraft and experiment modules. A flow diagram depicting the 
use and applications of T-013 flight experiment data is given as figure 14. 

CONCLUDING REMARKS 

This report summarizes the development of a manned space flight experiment to 
investigate the characteristics of crew-motion disturbances on a manned spacecraft. 

The data expected from this experiment will be used in the development of attitude con- 
trol systems for long-term manned space flight. An analytical crew -motion model, 
derived in the appendixes, may also be applied to the analysis of other man-motion 
phenomena and effects. Specific examples of potential applications include the use of 
selected limb movements to reorient a free-floating astronaut, the modeling of astronaut 
personnel vehicles (such as backpack maneuvering units), and work capability or task 
analysis (in conjunction with metabolic or energy expenditure). Verification of this math- 
ematical model through the Skylab Crew/Vehicle Disturbances Experiment should result 
in a valuable tool for research into a wide variety of manned space activities. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 23, 1971. 
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APPENDIX A 

SEGMENT MASSES, CENTERS OF MASS, AND LOCAL MOMENTS OF 
INERTIA FOR AN ANTHROPOMETRIC MODEL OF MAN 

By Charles T. Woolley 
Langley Research Center 

The derivation of a mathematical model of the human body, by which segment 
inertial properties and other biomechanical data can be computed is presented in this 
appendix. Figure 15 depicts a frontal view of the location of several principal skeletal 
anthropometric points used in the model derivation, along with locations of the body ref- 
erence planes. A lateral view of the remaining anthropometric point and the body ref- 
erence planes is shown in figure 16. 

An orthogonal coordinate system is used as a body reference system. This system 
is defined by the intersection of the three principal body planes shown in figures 15 and 16. 
The Xj^j-axis is formed by the intersection of the sagittal and transverse planes, with 
+Xmi directed toward the front of the body; the Yj^-j-axis is defined by the intersection 
of the frontal and transverse planes, with +Y]yu directed toward the right side of the 
body; the Z^jj-axis is defined by the intersection of the sagittal and frontal planes, with 
+Zmi directed toward the feet to complete the right-handed axis system. The origin of 
the coordinate system is located at the intersection of the three body reference planes 
as shown in figures 15 and 16. 

Figure 17 illustrates the relation of the body segments. Segment 1 is treated as a 
rigid member and results from the combination of three solids: the head is represented 
by an ellipsoid, and the upper and lower torso are represented by right elliptic cylinders. 

A lower arm segment (segments 4 and 5) is a combination of a sphere for the hand and 
a frustum of a right circular cone for the lower arm. Another compound segment is the 
lower leg (segments 8 and 9); this segment is a combination of two frusta of cones, one 
representing the lower leg, and another representing the foot. All other body parts are 
approximated by simple frusta of right circular cones. 

To use the mathematical model in crew-motion simulations, the relative locations 
of the segment centers of mass and the segment hinge points must be specified. Fig- 
ure 18 shows these locations, along with orientation of the coordinate axes. 

The model of the human body described in this appendix is a modification and exten- 
sion of a 15 -segment model developed by Hanavan (ref. 5). The current development 
reduces the number of segments to nine by the combination described above. It is recog- 
nized that the nine -segment model possesses possible limitations in use for crew- motion- 

26 


m ii i 


I 




Figure 15 . - Frontal view of htunan skeleton, showing location 
of body-reference planes and principal anthropometric 
points used in derivation of analytical model. 


Frontal Plane 



Figure 1 6.- Lateral view of human skeleton, showing 
location of one anthropometric point and the 
body reference planes. 
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Figure IT-- Nine-segment model of man. 


disturbance analyses. No account is taken of wrist, ankle, or head motions, nor of bend- 
ing of the spine. However, the hands and feet are relatively small masses when com- 
pared to other body segments, and would provide minimal disturbance effects if moved 
vigorously; whereas, the head and spinal motions involve larger masses which could 
affect T-013 experiment data. The T-013 motion sequences will be constrained to move- 
ments (insofar as possible) of only the nine principal segments. 

Twenty-four standard anthropometric measurements are used to dimension the 
model segments. Table V contains a list of these measurements, their abbreviations 
which will be used in the development to follow, and brief descriptions of the measure- 
ments. Complete descriptions of the measurements are found in reference 4. 


The man's total body mass must be distributed over each of the body segments. 
Barter (ref. 6) devised a set of regression equations which will be used for the distribu- 
tion process. These regression equations are given below: 


Mass of: 

Head, neck, and trunk = HNT = 0.47m^, + 5.44 kg 
Both upper arms s BUA = 0.08mJj, - 1.32 kg 
Both forearms s BFO = 0.04mlp - 0.23 kg 
Both hands s BH = O.Olm^ + 0.32 kg 
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TABLE V.- ANTHROPOMETRIC MEASUREMENTS 


Measurement 

Abbreviation 

Brief description 

Ankle circumference 

ANKC 

Minimum circumference of right ankle 

Axillary arm circumference 

AXLLC 

Horizontal circumference of upper arm at armpit 

Buttock depth 

BUTTD 

Front -to -rear dimension of body at buttocks 

Chest breadth 

CHESB 

Width of body (less arms) at level of nipples 

Chest depth 

CHESD 

Front -to-rear dimension of chest at level of nipples 

Elbow circumference 

ELBC 

Circumference of arm at elbow (arms extended) 

Fist circumference 

FISTC 

Circumference of fist measured over thumb and knuckles 

Forearm length 

FOARL 

Length of extended lower arm between elbow and wrist 

Foot length 

FOOTL 

Maximum length of foot from heel to toe 

Knee circumference 

GKNEC 

Circumference of right knee at kneecap (standing) 

Head circumference 

HE ADC 

Maximum circumference of head above eyebrows 

Hip breadth 

HIPB 

. Maximum width of body at hips 

Shoulder (acromial) height 

SHLDH 

Distance from floor to right shoulder (acromion), standing 

Sitting height 

SITH 

Height from seat to top of head, subject sitting erect with 
knees bent 90° 

Sphyrion height 

SPHYH 

Distance from floor to sphyrion (ankle), subject standing 

Stature 

STAT 

Distance from floor to top of head, subject standing erect 

Substernale height 

SUBH 

Distance from floor to substernale (lower breastbone), subject standing 

Thigh circumference 

THIHC 

Horizontal circumference of upper right thigh, standing 

Tibiale height 

TIBH 

Distance from floor to tibiale (knee), subject standing 

Trochanteric height 

TROCH 

Distance from floor to trochanterion (hip), subject standing 

Upper arm length 

UPARL 

Distance from acromion to radiale (shoulder to elbow), arm at side 

Waist breadth 

WAISB 

Minimum width of waist, standing 

Waist depth 

WAISD 

Front -to-rear dimension of body at waist 

Wrist circumference 

WRISC 

Minimum circumference of wrist, arm extended 
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APPENDIX A — Continued 
Segment Development 

The following section will develop the properties of the nine body segments. Each 
segment will be assigned a coordinate axis system, with its origin at the segment center 
of mass. Orientation of each coordinate system will be such that the segment X-, Y-, 
and Z-axis will be parallel to the Y MI -, and Z^-axis, respectively, when the body 

is in a normal standing position, with arms hanging vertically by the sides. The segment 
properties development follows that of Hanavan in reference 5. 

Hea d and torso (segment 1) .- As shown in figure 19, the head of the model man is 
represented by a right circular ellipsoid. The segment is oriented so as to produce a 



Figure 19. - Head of model man (right circular ellipsoid). 

circular cross section in a plane parallel to the transverse plane. The dimensions and 
properties of the head are given by 

R h = 0.5(STAT - SHLDH) 

RR H = HEADC/2ir 
SL h = STAT - SHLDH = 2R H 
Rjj = 0.5SLjj 
V H = 4irR H • (RR H ) 2 /3 
m H = (from Dempster, ref. 3) 

p H = m H /V H 

SI XX,H = °- 2m H [( R h) 2 + ( RR h) 2 ] 

SI YY,H = SI XX,H 

si zz,h - °- 4m H • (RRh) 2 

Right elliptical cylinders are used to simulate both the upper and lower torso. 
The orientation is such as to produce elliptical cross sections in planes parallel to the 
transverse plane. 
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The mass of the entire torso is found by subtracting the mass of the head from the 
mass of the head, neck, and trunk. A ratio of the densities is then used to find the mass 
of each of the upper and lower torsos. 

The upper torso is shown in figure 20. The dimensions and properties of the upper 
torso are given by 


where 


R ut = 0.5CHESB 
RR ut = 0.25(CHESD + WAISD) 
SL ut = SHLDH - SUBH 
SRut = 0.5SL^jij* 

Vut = ft ’ R^jt " rr ut * slut 


V LT = 77 ’ R LT ’ RR LT ' sl LT 

Put - ( HNT - “h)/( v ut + rV LT) 




Specific gravity^ 

r = 

Specific gravity 


~~ (from Dempster, ref. 3) 

u . y z 


(A2) 


m UT = P UT ' V UT 

SI XX,UT = m UT 3 ( r Ut) 2 + ( sl Ut) 

SI YY,UT = m UT 3 (R r ut) 2 + ( SL Ut) 


'12 


'12 


SI ZZ,UT “ m UT 


(r ut ) 2 + (RR UX ) : 


J 
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Figure 20.- Right elliptic cylinder representing the 
upper torso of model man. 
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Figure 21 shows the lower torso of the model. The dimensions and properties of 
the lower torso are given by 


R lt = 0.5HIPB ''v 

RR lt = 0.25(WAISD + BUTTD) 

SL lt = SITH - (STAT - SUBH) 

SRjj'p = 0.5 SL lt 


v lt = Jr • r lt • rRl T • sl lt 

= HNT - rtipj - m-yrp 
P LT = m LT/ V LT 


XX, LT ~ LT 


{3(R LT ) 2 +(SL LT ) 2 j/l2 

SI YY,LT = m LT [3(RR LT ) 2 + ( SL Lt ) 2 ]/ '*■ 

__ / 

SI ZZ,LT = m LT^( R LT) 2 + ( RR Lt) 2 






Figiire 21- - Right elliptic cylinder representing the 
lower torso of model man. 


(A3) 


These first three solids are combined to form the first segment of the model man. 
The resultant segment is shown in figure 22. The dimensions and properties of the total 
torso segment are given by 
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DELSH = SITH - STAT - TROCH^I 


SLp = SLjj + SLyp + SL^rp 
SMLp = SLp - DELSH 


m t = mjj + m TT rp + m TT = HNT 


UT 


LT 


Vi = V H + V UT + V LT 
Pl = mp/Vp 


(A4) 



Figure 22.- Segment 1 of model man (from combination of 
head, upper torso, and lower torso segments) . 


The distance from the top of the head to the center of mass for this multibody sys- 
tem is 


m 


SMRp = 


H ' % 


+ m 


UT 


(SL 


H 


+ SR 


UT, 


+ m^T (SLjj + SLurp + SR j 


HNT 


(A 5) 


The total local moment of inertia about each of the segment reference axes is 
given by the parallel -axis theorem 
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SI XX,1 = SI XX,H + SI XX,UT + SI XX,LT + m H( SR H " SMRj ) 2 

+ ni-yrp ^SLjj + SR-yrp — SMR-j^ + mj^y ^SLj^ ^ R LX ~ SMR^ (A6) 

SI YY,1 = SI YY,H + SI YY,UT + SI YY,LT + m H( SR H " SMR l ) 2 

p p 

+ m xJT(S^- , H ®^UT " SMRjJ + SL^jrp "* (A7) 

SI ZZ,1 = ZI ZZ,H + SI ZZ,UT + SI ZZ,LT ( A8 ) 

Conical segment properties .- In all following development, the properties of conical 
segments will apply. To simplify the discussion, all necessary conical segment proper- 
ties are listed. 


V = 

SL • tt(r 2 + RR 2 + R • RR)/3 

~\ 

M = 

RR/R 




a = 

1 + fi + 

M 2 



V = 

(l + 2 m 

+ 3ju 2 )/ Ict 



AA 

= 9(l + 

9 q 4 

fl + /jA + JUl + fl^ 

)/ 20 ttct 2 

BB 

= 3(l + 

4 JU + 10p2 + 4^3 

+ m 4 ) ; 

/80ct 2 

SR 

= rj • SL 




SIxx - AA 

• m 2 /p • SL + BB 

* m • 

SL 2 


Slyy = SIxx 


SI^z = 2AA * m 2 / p • SL 


(A9) 




Upper arm (segments 2 and 3) .- The upper arm segment is represented by the 
frustum of a right circular cone. This segment is represented in figure 23. The prop- 
erties and dimensions of the upper arm are given by 
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For this segment, the local moments of inertia and remaining properties are those 
of a conical segment. 



Figure 23-- Upper arm of model man (segments 2 and 3 )- 

Lower arm and hand (segment s 4 and 5) .- The hand of the model is represented by 
the sphere shown in figure 24. The properties and dimensions of the hand are given by 



Figure 24 . - Hand of model man . 
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r HA = FISTC/2 tt 
rr HA = r HA 
sl ha = 2R HA 
^HA = 0,5 
m HA = 

sr ha = r ha 

P HA = 3m HA/ 4?r ' R HA 3 
SI XX,HA = °- 4m HA ' R HA 2 
SI YY,HA = SI XX,HA 
SI ZZ,HA = SI XX,HA 


(All) 


The lower arm of the model is represented by a frustum of a right circular cone. 
The orientation is such as to produce a circular cross section in a plane parallel to the 
transverse plane. The segment is shown in figure 25. The dimensions and properties 
of the lower arm are given by 



Fi gur e 25 . - Lower arm of model man . 


R la = ELBC/27T 
RR t a = WRISC/2tt 

> 

sl la = FOARL 

m^ A = 0.5BFO 


(A12) 


Local moments of inertia and remaining properties are those of a conical segment. 


37 



I I ■ 1 1 I 


APPENDIX A — Continued 

In figure 26, the assembled lower arm and hand segment is shown. The hand seg- 
ment and lower arm segment are combined as the total lower arm segment used in the 
nine -piece model. The dimensions and properties of the combined segment are given by 



SL 4 - SLjja + SL la 

m 4 = m HA + m LA 
v 4 - Vha + V LA 
P 4 = m 4 /V 4 


(A13) 


Figure 26.- Segments 4 and 5 of model man 
(from combination of hand and lower 
arm segments). 


The distance from the hinge point to the center of mass for this multibody system 
becomes 


SMR4 = 


m HA( SL LA + SR Ha) + m LA ' SR LA 

m4 


(A14) 


The total local moment of inertia about each of the segment reference axes is 
derived by the parallel-axis theorem, and yields 

SI XX,4 ~ SI XX,HA + SI XX,LA + m HA( SL LA + SR HA " SMR 4) 2 + m LA( SR LA “ SMR 4 ) 2 

(A 15) 

SI YY,4 = SI XX,4 (A16) 

SI ZZ,4 = SI ZZ,HA + SI ZZ,LA (A17) 

Upper leg (segments 6 and 7).- The upper leg of the model is represented by the 
frustum of a right circular cone. Orientation is such as to produce a circular cross 
section in a plane parallel to the transverse plane. The segment is illustrated in fig- 
ure 27. The properties and dimensions of the upper leg are given by 
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R 6 = TfflHC/277 

RR 6 = GKNEC/2tt 

SLg = STAT - SITH - TIBH 

SMLg = SLg + DELSH 

SMRg = SRg + DELSH 

mg = 0.5BUL 


(A18) 


Figure 27.- Upper leg of model man 
(segments 6 and 7). 

Local moments of inertia and remaining properties are those of a conical segment. 

Lower leg and foot (s egments 8 and 9) .- The lower leg of the model is represented 
by the frustum of a right circular cone. The orientation is such as to produce a circular 
cross section in a plane parallel to the transverse plane. Figure 28 depicts this subseg- 
ment. The properties and dimensions of the lower leg are given by 



Figure 28.- Lower leg of model man. 


R ll = GKNEC/2tt 
RRLL = ANKC/2tt 

> 

SL ll = TIBH - SPHYH 
m LLj = 0.5BLL 


(A19) 
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The local moments of inertia and remaining properties are those of a conical 
segment. 

The foot of the model is also represented by the frustum of a right circular cone. 
The orientation is such as to produce a circular cross section in a plane parallel to the 
frontal plane and a trapezoidal cross section in planes parallel to the transverse and 
sagittal planes (fig. 29). The dimensions and properties of the foot segment are given by 

R f = 0.5SPHYH ^ 

RR F = 0.403R F (from Dempster, ref. 3) 

SL f = FOOTL \ (A20) 

m F = 0.5BF 

77p = 0.429 (from Dempster, ref. 3) 

J 

The local moments of inertia and the remaining properties are those of a conical 
segment. 



The total lower leg segment of the nine -piece model is composed of the two pre- 
ceding segments, as shown in figure 30. The dimensions and properties of the combined 
lower leg segment are given by 


SL g = SL ll + SPHYtP 

V 8 = v LL + V F 

m 8 = m LL + m F ( 


(A21) 


PQ = m 8/ V 8 
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Figure 30.- Segments 8 and 9 of model man (from combination 
of lower leg and foot segments). 


The distance from the hinge point to the center of mass for this multibody system 
is given by 


SMRg = 


m LL ' SR LL + m F( SL LL + %) 
m 8 


(A22) 


The total local moments of inertia about each of the segment reference axes follows 
from the parallel-axis theorem 


®^XX,8 = SI XX,LL + SI XX,F + m LL( SR LL ' SMR s) 2 + m F( SL LL + R F " SMR s) (A23) 

Slyy g = Slyy^LL + Slyy p + m LL (S^ll - SMRg)^ + mp^SLLL + Rp - SMRg)^ (A24) 


SI ZZ,8 " SI ZZ,LL + SI ZZ,F 


(A25) 


Comparisons and Extensions 

The equations defining segment properties have been incorporated into a computer 
program, called MAID (for Man’s Analytical Inertia Determination), which computes total 
body inertias and center-of-mass location from the individual segment properties. 

Inputs required are the subject's mass and the 24 anthropometric measurements listed 
in table V. Comparison of MACD-computed inertias and center-of-mass location with 
experimentally measured inertias and center-of-mass locations was performed, using 
data from Santschi, DuBois, and Omoto (ref. 2). Anthropometric data from 66 Air Force 
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flight personnel, contained in reference 2, were used as the input to program MAID. The 
computed data from MAID were then compared with the measured inertial and center -of - 
mass properties of the same 66 subjects, with the results presented in table VI. 

TABLE VI.- COMPARISON OF MEASURED AND COMPUTED 
INERTIAL PROPERTIES FOR ENTIRE MAN 



Santschi et al. (ref. 2) 
(measured) 

MAID 

(computed) 

Average error, 
percent 


High 

Low 

Average 

High 

Low 

Average 

(ref. 2) 

Local moment of 
inertia about: 








X-axis, 

kg-m^ 

18.63 

8.52 

13.00 

18.51 

8.11 

13.05 

-0.5 

Y -axis , 
kg- m 2 

16.50 

7,49 

11.58 

17.23 

7.61 

12.30 

-6.2 

Z-axis, 

kg-m 2 

1.97 

.77 

1.28 

1.81 

.72 

1.10 

+13.8 

Location of c.m. 
along Z^-axis, 
m 

0.072 

0.059 

0.066 

0.071 

0.060 

0.066 

0.0 


Over the range of the 66 test subjects, there appears to be good agreement between 
the MAID and experimental results for X- and Y-axis total inertia values and for center- 
of-mass locations. The larger difference in the Z-axis inertia comparison has been dis- 
cussed by Hanavan (ref. 5), who performed similar comparisons with his model. Since 
the moments of inertia for the Z-axis are an order of magnitude smaller than those for 
the X- and Y-axis, they are subject to larger relative experimental error. It was also 
pointed out by Hanavan that, upon inspection of documentary photographs taken during the 
measurements reported in reference 2, it appeared that subjects were not held precisely 
in position. 

The average Z-axis inertia calculated in reference 2 is quite large when compared 
to that computed by MAID. Since any misalinement with respect to the Z-axis by the sub- 
ject would cause the measured inertia to be larger than the actual value, the explanation 
presented by Hanavan appears to be reasonable; the analytical model is not subject to 
these misalinement errors. The large deviation from experiment for the Z-axis inertial 
properties should therefore not invalidate the analytical model. 

Reference 2 noted that there appears to be a very high degree of correlation between 
body inertia and total mass. It seemed reasonable to assume that individual segment iner- 
tial properties could also be correlated with body mass. Based on the work by Santschi, 
DuBoise, and Omoto, a linear least-squares analysis was performed to determine 
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regression equations for local moment of inertia versus total body mass for each seg- 
ment. The data from the 66 flight personnel, reported in reference 2, were used in the 
least-squares fit analysis. Results of the least-squares-curve fit are presented as fig- 
ures 31 to 35. These figures present the experimental data points calculated from the 
subject dimensions given in reference 2, computed data points for five different percentile 
cases, and the lines of regression (solid lines), their standard deviations (broken lines), 
and algebraic equations. Segment moment of inertia is plotted against total body mass 
for each figure. 

The five percentile cases, denoted by inverted triangles on the figures, represent 
composite, and not actual, subjects. For example, 5th-percentile values were taken from 
each set of the 24 anthropometric measurements in order to make up the 5th-percentile 
man. The 5th-, 25th-, 50th-, 75th-, and 95th-percentile cases were computed for each 
segment inertia and plotted against corresponding subject mass percentile values (60.2 kg, 
67.5 kg, 73.6 kg, 80.35 kg, and 91.2 kg, respectively). It can be seen from figures 31 to 35 
that the linear-curve fits were very good for all segments as indicated by the small stand- 
ard deviations. Also, the agreement of the five percentile cases with the regression equa- 
tions was close. It is thus apparent that the regression equations can be a useful tool in 
computing segment inertial properties, with only a knowledge of total body mass of a par- 
ticular subject. 
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Figure 31- - Variation of segment moment of inertia with man's 
mass for segment 1 of model man. 
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Figure 31- - Concluded. 
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• Subject data (from ref. 2) 
V Percentile eases (computed) 
T.ine of regression 


X ZZ * 0.0000335 m - 0.001041 
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Figure 33 . - Variation of segment moment of inertia with man*s mass 
for segments 4 and 5 of niodel man. 
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Figure 34.- Variation of segment moment of inertia with man’s mass 
for segments 6 and 7 of model man. 
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• Subject data (from ref. 2) 
V Percentile cases (computed) 
Line of regression 
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Figure 35 -- Variation of segment moment of inertia with man's mass 
for segments 8 and 9 of model man. 
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DEVELOPMENT OF A DYNAMIC ANALYTICAL MODEL OF MAN 
ONBOARD A MANNED SPACECRAFT 

By Peter R. Kurzhals and Robert B. Reynolds 
Langley Research Center 

This appendix utilizes the model man developed in appendix A, and derives applied 
force and torque equations corresponding to motions of this model. These torque and 
force equations are incorporated in the rigid-body equations of motion for arbitrary 
spacecraft. Both the disturbance torque due to crew activity and the effects of these 
torques on the spacecraft attitude and rates are determined by specifying the motion of 
the mathematical model within the spacecraft. 

Reference Coordinates 

The reference coordinates for the dynamic-model development are illustrated in 
figure 36. The spacecraft angular motion is defined by a set of X,Y,Z of body-fixed 
axes which rotate with respect to a set Xj,Yj,Zj of intermediate axes. The origin of the 
X,Y,Z system is selected as the spacecraft center of mass, without any crew (or model 
man). 

The intermediate axes are free to translate without rotation in inertial space, and 
always remain parallel to a set Xp,Yp,Zp of inertially fixed axes. The attitude of the 
spacecraft is described by Euler angles which relate the X,Y,Z- and Xj,Yj,Zj-axes. 

The rotation of the model man is specified by two sets of axes whose origins are 
coincident with the model center of mass. The Xjyj,Yjyj,Zjyj-set is parallel to the principal 
axes of segment 1 of the model (torso and head), and rotate relative to the x mS’ Y MS’ Z MS" 
axes, which remain parallel to the X,Y,Z-axes. The +X M -axis points toward the front of 
the model; the +Yj^-axis points toward the model's right side; and the +Z]yj-axis points 
toward the model's feet. 

The motion of the model man’s center of mass is related to the man's lower pivot 
points (between torso and upper legs) by rectangular coordinates in a set 
of body-fixed axes. These axes remain parallel to the X M ,Yjyi,Z]y[-axes and have their 
origin at a point halfway between the hinge points noted in figure 18. The total translation 
of the model's center of mass is given by the sum of the translation of the X MI >YmI’ Z MI - 
origin and the translation of the model’s mass center with respect to this origin. Mass 
center coordinates must be expressed as rectangular coordinates in the X,Y,Z-axis 
system. 
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Figure 36.- Reference coordinates used in the development of 
the dynamic model of man. 

The angular position of segment 1 is described by three Euler angles, shown in fig- 
ure 37, which relate the x M> Y M’ Z M -axes to the x MS’ Y MS> z MS -axes anc * hence to the 
X,Y,Z -space craft -fixed axes. The angular motion of all other body segments is described 
by two Euler angles per segment, which relate the moving segments to axes parallel to the 
x M> Y M’ Z M" axes of se g ment 1, as shown in figure 38. 

The model’s motion with respect to the spacecraft is defined by specifying the time 
variation of the origin of the Xj^Yj^jZ^ - system (or, eventually, the model center of 
mass) and of the Euler angles for segments 1 to 9. 

Spacecraft Euler Angles 

The inertial spacecraft attitude is defined by three modified Euler angles, illus- 
trated in figure 39, which determine the relative motion between the X,Y,Z- and Xj,Yj,Zj- 
axes. These modified Euler angles are designated as if/, 9, and cp, and result from 


52 




APPENDIX B - Continued 



Figure 39*“ Euler angles used to describe spacecraft rotations. 
Rotation sequence is cp. 


t tt 

three consecutive rotations about the Zj,YpXj -axes. From figure 39, the time deriva- 
tives of the Euler angles become 


cp- O x + Qy tan 0 sin cp + £ 2 Z tan 0 cos cp ^ 
0 = £ 2y cos cp - £l z sin cp ’ 

xfy = £2 Z cos cp sec 0 + fiy sin cp sec 0 


(Bl) 


The Euler angles can be found by numerical integration of equations (Bl), after the body 
rates £2y, and £2 Z have been derived by the numerical integration of the equations 

of motion for the spacecraft. Time histories of both the Euler angles and the body rates 
are used to describe the angular motion of the spacecraft. 


Equations of Motion 

To develop the equations of motion for the spacecraft and man, the following nota- 
tion is introduced. Any 3X1 vector, such as the spacecraft angular velocity <^T}> and 
the man’s center of mass position <(Y}>, is designated by braces A. Any 3X3 matrix, 
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such as the matrix |jg] of the spacecraft inertia tensor, is designated by brackets 
A vector symbol which is enclosed by brackets instead of braces represents the skew 
symmetric matrix corresponding to that vector. For the spacecraft angular velocity 




(B2) 


this skew symmetric matrix is 



-O z 

0 

Ox 


Oy 

-Ox 

0 


A similar matrix can be constructed for any vector. 


(B3) 


The moment equation in the geometric spacecraft coordinate axis system X,Y,Z 
may now be written as 


{M E } + {M M } = + ^L[r]{F s } 


(B4) 


where the external moment , the moment <(M]y0> exerted by man, and the radius 

vector {yy to the man's center of mass are all referred to the spacecraft geometric 
origin. The spacecraft angular momentum ^Hg^> and the total force <^Fg^> are referred 
to the changing composite mass center of the spacecraft and the moving man. 

The moment exerted by crew activity is derived by considering the man as a sepa- 
rate dynamic system whose motion is completely specified in spacecraft coordinates. 

This leads to a moment expression of the form 


( m m) = 


d ( H M) 


dt 


m M [ r ]- 


dt 2 


(B5) 


where is the angular momentum of the man referred to his mass center and {jEl} 

is his mass center radius in fixed (inertial) coordinates. The second term on the right- 
hand side of equation (B5) represents the torque when the man is taken as a point mass 
nijyj located by the radius vector in spacecraft geometric coordinates. The quan- 

d 2 /R} 

tity -mjyr — represents the force that is applied to the spacecraft at the point where 

dt 2 

the man is in contact with it. Components of this force are measured for the subject of 
experiment T-013 by the force measuring units. 
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The acceleration of the man’s mass center is given by 

d2 fo) _ /p T d2 ( f ) 

dt 2 dt 2 

= {Rq} + & + 2 M (r) + [n]{r> + [«] [o] {r} (B6) 

where <(R 0 ^> * s P° s Ri° n vector in fixed coordinates to the spacecraft geometric 
origin. Since 


<*S> - <Je} - “s({Rq> + ^-^r) 


(B7) 


where denotes any external force, the torque produced by the total spacecraft 

force aR d the torque produced by man taken as a point mass can now be combined 


as 


1 — | ip* \ ^{*5} /*• \ m M ^ © 1 r ^ 

-mH K r o> + - (Ro> - -st -££ + ife <?e} 


2 < 

S dt" 

2 j 

dt 2 




(B8) 


Substitution of equations (B5) and (B8) into equation (B4) yields 

r M -> d ( H M> _ d { H S> m M( m S - m M) r ] d2 ® m Mr rl r , 

{. m E> - — ^ dt~ + me L r J^2~ + IEtI/JV 


f e} 


Let 


and 


to get 


T = - 


Q = 


1 ( h m} 

dt 

^(mg - m M ) 

m c 


( m eM t > dt 

{h s } , [i s ] (n) 




m 


M 


dt 2 + m S 


H ( f e> 


(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 
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and 


M- 


*xs 

-I XYS 

_I XZS 

-I XYS 

X YS 

-I YZS 

_I XZS 

_I YZS 

*ZS 


(B14) 


In the expanded form the equations of motion represented by equation (B12) become 

{Me} + {T} = [I S ] {§ + [O] [Is] {n} + Q[r] (0) + 2 [o] {r} 

For this development, it is convenient to introduce a total angular momentum 
given by 

,4} 


(B15) 


(H X } = {H S } + Q[r]. 


dt 


(B16) 


which is equal to the angular momentum of the spacecraft and the moving man taken as a 
point mass. Equation (B16) can be rewritten in terms of the familiar "rigid-body" angu- 
lar momentum, denoted by which is formulated in terms of the instantaneous 

inertia of the spacecraft and point -mass man and the spacecraft angular rates. The 
result is 


where 


(H t ) = {H RB } + Q[r]{r} 

{%B> = [%] 0) + qHM 0} = [i] ® 


with 


M- 


-I 


XY 


-I 


XZ 


-I 


YZ 


I X 

_I XY r Y 
_I XZ _I YZ *Z 

Substitution of equations (B16), (B17), and (B18) into equation (B12) yields 


(B17) 

(B18) 

(B19) 


{m e } ♦ {t} = + ^[r] {Fe} 

+ ^H{Fe} (B20, 
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or 


{m e } + {t} = [i]{n} + [i]{n} + [n][i]{n> 

+ Q([r](r} + [ n ][r](0) ♦ {F E } (B21) 

as the final equations of motion. 

Spacecraft Inertia and Mass Center Coordinate Terms 

The spacecraft inertia terms in equations (B19) and (B21) may be written in 
terms of the spacecraft inertias (excluding the model man), the mass factor Q, and 
the Cartesian components of the man's center of mass radius vector as 


J x = I xs + Q(y 2 


+ z 




Iy = X YS + Q(x 2 + z 2 ) 

iz = izs + q ( x2 + y 2 ) 

I XY = I XYS + Q ( x y) 

^XZ = *YZS + Q( xz ) 

- *yzs + Q(y z ) 

The time derivatives of the inertia terms, which appear in equation (B21), become 
I x = 2Q(yy + zz)^ 
fy = 2Q(xx + zz) 


(B22) 




(B23) 


i z = 2Q(xx + yy) 

ixy = Q( x y + yx) 

i xz = Q(xz + zx) 

iyz = Q(y z + zy) 

-X 

The radius vector to the center of mass of the mathematical model of man is 

specified by 

{r} = {r T } + (r cm ) (B24) 

where ^Tr0, is the radius from the spacecraft geometric origin to the origin of the 

X MP Y MI’ Z MT axis s y stem > and ( r cm) is the radius vector from the 

origin to the mass center of the model man. Remembering that <{V}> must be expressed 

in X,Y,Z- coordinates, |r cm ^> is expressed in that system by the transformation 
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f x ]Vn| 

(pi) \ ym / = [ D l] <C r Ml} 

l z MU 


(B25) 


The Euler angle transformation matrix [pi] for segment 1 of the mass center 
coordinates x ffl , y M£ , and Zj^j in the X MI> Y MI’ Z Mr system are defined in subse- 
quent sections of this appendix. 

Time derivatives of {V} are needed in the equations of motion and are given by 

O'} = (Ft) + M {5 Ml} + [ D l]OMl} < B26 ) 

and 

0} = (r T ) ♦ [&l]pm)+ + [D^V) (B27> 


Segment Transformation Matrices 

Quantities measured relative to the x m> y M>Zm" or x MI> Y MI»ZMI~ axes ma y t> e 
referred to the X,Y,Z -system by means of the transformation matrix (p{J; likewise, 
quantities measured relative to the segment coordinate axes Xj,Yj,Zj (j = 2, 3, 4, . . .,9) 
may be referred to the x M> Y M> z M -axes by means of the transformation matrices [DjJ , 
where j = 2, 3, . . .,9. 

Matrices required for the transformation process are (j = 1, 2, . . 9 in each case) 



i d 12 

jd 22 

3 d 32 


3 d 13 

jd 23 


i d 33 


(B28a) 


jdn 

3 d 12 

3 d 13 

J d 2i 

3 d 22 

3 d 23 

jd 31 

3 d 32 

3 d 33 


(B28b) 



3 d 12 

jd 22 


3 d 13 

jd 23 


3 d 31 j d 32 3 d 33J 


(B28c) 


59 



APPENDIX B - Continued 


jd 


-1 

11 


H- 


id 21 

IjdgJ 


jd 12 

Jd 22 

j d 32 


^ d 13 

id 23 

jd 33 


For segment 1 


N- 


cos yjg cos y^ 2 

cos y^g sin y^ sin y^ 

cos sin y^ 2 cos ^ 


- sin y^g cos y^ 

+ sin y^g sin y^ 

sin cos y ^ 

sin y^g sin y ^ sin y^j 

sin y^g sin y^ cos y^j 


+ cos y^g COS y^j 

- cos y^g sin y^ 

-sin y 12 

cos y^2 sin y^ 

cos yjg cos y^i 


[*■]- 


For segments 2 to 9 

'cos y. 2 sin y- 2 sin y.j sin y_ 2 cos y.j' 
0 cos y^i -sin y^ 

-sin y -2 cos y^ 2 sin y^ cos y^ cos y^ 


(j = 2, 3, . . 9) 


For all segments 






d ( Jd n) 

d (i d i 2 ) 

d 0 d i 3 ) 


dt 

dt 

dt 


d (i d 2 i) 

d 0 d 22 ) 

d 0 d 23 ) 


dt 

dt 

dt 


d 0 d 3 l) 

d 0 d 32 ) 

d 0 d 33 ) 


L dt 

dt 

dt J 


d 2 (i d n) 

d 2 (i d l 2 ) 

d 2 0 d 13 ) 

dt 2 

dt 2 

dt 2 


d 2 (j d 2 i) 

d 2 (j d 22 ) 

d 2 0 d 23 ) 

dt 2 

dt 2 

dt 2 


d 2 (i d 3 i) 

d 2 0 d 32 ) 

d 2 0 d 33 ) 

dt 2 

dt 2 

dt 2 

J 


(j = 1, 2, . . 9) 


(j = 1 , 2 , . . 9 ) 


(B28d) 


(B29) 


(B30) 


(B31) 


(B32) 
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[ D i] _1 * [ d j] 


(j — 1> 2, . . .,9) 


(B33) 


It should be noted that |~DjJ is needed only for j = 1 (i.e., for segment 1). 

Model and Segment Mass Centers in Model Coordinates 

By assigning center of mass and pivot -point coordinates to the mathematical model 
of figure 17, one arrives at the system shown in figure 18. 

The segment center -of -mass coordinates, in the -system can be read- 

ily derived from figures 18 to 30 with the result 


1MI 

y lMI 

Z 1MI 

X 2MI 

y 2MI 

Z 2MI 


r 


mmi 


y 4MI 


4MI 


■ 


0 

0 




(B34) 


-SMLj + SMRj 


- 


X 3MI 

{ y 3MI ) = { 
Z 3MI 


r 


SMR2 sin y 22 cos 7 21 
Rut + R2 - SMR2 sin 72^ 

-SMLj + S LI j + SMR 2 cos 722 cos y 21 

SMR2 sin 2 cos r 31 
“RuT - R 2 “ SMR2 sin 7 2^ 
-SMLj + S^u + SMR 2 cos ygg cos y^^ 

SMLo 




(B35) 




(B36) 


"A 


= 


x 2m + SMR 4 Sin y 42 cos y 41 
SML 2 

r UT + r 2 + s&m^( y 2MI " R UT " R 2y " SMR 4 sin 7 41 

/ sml 2 \ sml 2 

(-SMLi + SL H )(l - z 2MI + SMR 4 cos y 42 sin y 41 




(B37) 
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c ^ 

X 5MI 


/ 

SMLo 

__ x 3MI + SMR 4 sin y 52 cos y &1 

y 5MI 

>=< 

SMLo . 

-R ut - R 2 + gMR2 (y3MI " R UT " R 2) - smr 4 sin r 51 

Z 5MI 
V- J 


. ./ SMLo\ SMLo 

(" SML 1 + SL h)^ - SMR 2 j + SMR 2 Z 3MI + SMR 4 cos r 52 sin y 51 


(B38) 


C ^ 

X 6MI 


r ^ 

SMRg sin yg 2 cos yg^ 

\ y 6MI 

i 

)-< 

r LT - r 6 “ SMRg sin yg-A 

Z 6MI 

V y 


SMRg cos y g2 cos y 61 


(B39) 


f x 7Mll 


r 

SMRg sin yrj 2 cos 




< y 7M i ) = \ - r lt + r 6 “ SMR 6 sin r 71 > 


(^MlJ 




SMRg cos y 72 cos y ?1 J 


(B40) 


r 

X 8MI 


SMLg 

SMRg 


X 6MI + SMR 8 sin ^82 


cos y 81 







SML fi t 

R LT ■ R 6 + §^R^( y 6MI ' R LT + R 6j - SMR 8 sin >W 


Z 8MI 

V. J 




SMLg 

SMRg 


Z 6MI + SMR 8 cos ^82 


cos y 81 


(B41) 
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r 

x. 




9MI 


y 9MI 


r 


SML 


A 


SM^ X 7MI + SMR 8 sin r 92 cos y gi 


SML 


> = < " r LT + R 6 + si^( y 7MI + r lt " R e) “ SMR 8 sin Y 91 > 


9MI 
v. j 


(B42) 


SML 




Z 7MI + SMR 8 cos ^92 cos ^91 


J 


Time derivatives of these coordinates can be easily found remembering that segment 
dimensions are constants for a given model man. 

The required position coordinates (denoted by a bar) of the overall mass center are 
expressed as 


MI 


9 

z 

= j=l 


9 


m j X jMI 


m r 


y MI 


Z “i y 3 

= i=i 


9 


]MI 


Z m i : 


m n 


J MI 




'jMI 


m r 


(B43) 


are 


where m^, total mass of the man, is nij. Time derivatives for equations (B43) 

3=1 

easily found, with the values of m T and mj being constant. 

The mass center coordinates ( x M ,y M ,z M ) for the jth segment (j = 1, 2, . . .,9) are 
found by 


r "\ 

x 


‘Mj 

y MJ 
z 


"\ 


x jMI " X MI 


( ~ \ y jMI " y MI 


J Mj 


J 




z jMI Z MI 


or 


{ r Mj} = { r jMl} - {5 Ml} 
Derivatives are 

'O’Mj} = {i’jMl} - {^Ml} 


and 


{*M]) = {r jM l} - Qm} 


(B44) 


(B45) 

(B46) 

(B47) 
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Specification of functions for the Euler angles yj and the translation of the 
^MD^MD z MI _or i£i n are re( l u i re d to express the motion of the model man. 


Disturbance Torques for Man’s Angular Motion 

The disturbance torque resulting from angular motions of the model man can be 
readily derived in the Xj^,Y]y[,Zjyj-axes system, which moves with the man's torso and 
has its origin at the man's center of mass. The torque in the spacecraft coordi- 
nate system is related to the torque acting on the man in the XjyjjY^jZ^-system 

by the transformation 

{t} = -[d 1 ]{Tm} (B48) 

where 

d ( H M} 


( T M> 


dt 


(B49) 


and is now expressed in the man's coordinate system. The man's angular momen- 

tum can be written as 


9 9 

( H M> = [%]{%[} + ^ [ D j] [Ij] [ D j] 1 {°j} + ^ [ m ][ r j] { r j} (B50) 

j=2 j=l 

where is the instantaneous inertia matrix of the man, considering all segments, 

pj] are the principal inertia matrices for the limbs moving relative to the torso, 
is the total angular velocity vector of the torso, and are the angular velocities of 

the body segments 2 to 9 relative to the torso. 

Substitution of equation (B50) into equation (B49) gives 





9 


+ \ m 3([ r i]^j) + [° M ] [ r j] C r J>) 


(B51) 


Expressions for the angular velocity vectors for the man's limbs and torso are needed 
for the evaluation of equation (B51). The total angular velocity of the torso is given by 
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r 


y n - y 13 sin y^ 

{Pm} = [ D l] " 1 {°} + <j Y 12 cos rn + r 13 sin r u cos y 12 

y 13 cos y 12 cos y n - y 12 sin yjj 

while the angular velocities of the limbs relative to the torso are 

Yu cos y. ^ 




(B52) 


fi-i — 


j “< 


]1 r ]2 
^]2 




(i = 2, 3, . . 9) (B53) 


-Yu sin y 


]2 


The moments and products of inertia for the model man are given by the matrix 
relation 

9 


where 


[I M ] - (pi] * [ill]) + J, ([Djpjpj]' 1 + pi]) 

3-2 


[i M ] - 


(B54) 


M- 


-I 

x xj 

0 

0 


[ XM 

-I YXM 

-I ZXM 

XYM 

t ym 

-I ZYM 

XZM 

-I YZM 

J ZM 


(B55) 


r Yj 


0 

0 


*Zj 


(i = 1, 2, 


9) (B56) 


and 




2 2 
y Mj + z Mj 

-x Mj y Mj 

" x Mj Z Mj 


= m i 

~ x Mj z Mj 

2 2 
x Mj + z Mj 

" y Mj z Mj 



" X Mj z Mj 

“ y Mj z Mj 

2 2 
X Mj + y Mj 


(j = 1, 2, 


9) (B57) 


Values of |jbjj and XMi’ y Mj’ z Mj have been developed previously. 

Use of the model inertias from equation (B54) and their time derivatives, along with 
the angular velocities from equations (B52) and (B53) and their time derivatives allows 
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the determination of {T m J, in equation (B51), and hence the torque {¥} exerted on the 
spacecraft by man in equation (B48). The total moment applied to the spacecraft can then 
be computed by using equation (B5). 
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